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ABSTRACT 


Tiie  clitnatoioey  of  clear-air  turbulence  is  defined  herein  as  the 
likelihood  that  an  aircraft  or  missile  will  encounter  turbulent  air  at 
a  given  locality,  altitude,  and  time  of  year.  Turbulence  data  of  three 
types  were  used  in  this  study;  these  include  obaervations  by  instrumented 
research  aircraft,  balloon  tracks  measured  by  FPS-16  radar,  and  turbu¬ 
lence  reports  made  by  pilots.  The  measurements  that  have  been  made  by 
research  aircraft  show  good  relationships  between  turbulence  and  certain 
aspects  of  mesoscale  atmospheric  structure,  but  the  data  are  too  limited 
in  number  to  permit  broad  generalizations.  The  FPS-16  tracks  of  rising 
Jimsphere  and  Rose  balloons  were  obtained  during  studies  of  detailed 
wind  profiles.  We  investigated  their  potential  value  in  identifying 
turbulent  layers.  In  general,  the  existing  data  are  too  noisy  to  serve 
this  purpose;  however,  further  special  trials  are  recommended.  The  sub¬ 
jective  turbulence  reports  from  pilots  collected  during  special  five-day 
reporting  periods  <^omprise  by  far  the  largest  volume  of  data  available. 

Meteorological  conditions  for  these  periods  were  analyzed  by  computer 
from  standard  rawinsonde  data  and  were  correlated  with  the  turbulence 
reports.  The  correlations  show  that  the  vertical  vector  wind  shear 
corresponds  most  closely  to  turbulence  frequency  determined  from  the 
pilot  reports.  Little  additional  reduction  of  variance  in  the  turbulence 
frequencies  is  achieved  by  including  other  meteorological  factors.  Op¬ 
timum  multiple  regression  equations  between  turbulence  frequency  and  the 

mean  and  standard  deviation  of  the  vertical  vector  wind  shear  were  ob-  • 

tained.  In  summer  a  different  regression  equation  was  found  than  in 
other  seasons. 

These  turbulence  observations  taken  in  toto  are  too  few  to  permit 
a  direct  computation  of  turbulence  frequency;  therefore  an  indirect 
method  was  used  to  obtain  the  turbulence  climatology.  This  involved 
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applying  the  regression  equations  to  existing  statistics  of  wind  shear 
over  the  United  States.  Tnese  wind-shear  compilations,  prepared  almost 
ten  years  ago,  appear  tc  rf  questionable  reliability  in  the  upper 
troposphere  and  stratosphere  due  to  a  large  proportion  of  missing  wind 
observations  under  conditions  of  high  wind  speeds  aloft.  Due  to  uncer¬ 
tainties  in  the  regression  equations  and  shear  statistics,  the  deduced 
turbulence  climatology  (given  by  seasons  for  levels  between  approxi¬ 
mately  20,000  and  45,000  feet)  must  be  considered  a  first  estimate. 

Recommendations  are  made  that  an  up-to-date  wind-shear  climatology 
be  computed  for  the  United  States,  and  that  consideration  be  given  to 
developing  a  balloon-borne  turbulence  sensor  to  augment  data  from 
research  aircraft  and  airline  pilots. 
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I  INTRODUCTION 


The  purpose  of  this  study  was  to  develop  methods  for  estimating  the 
probability  that  an  aircraft  will  encounter  turbulent  air  at  a  given 
locality,  altitude,  and  time  of  year.  The  discussion  pertains  to  clear- 
air  turbulence,  which  is  perhaps  more  accurately  described  as  non-storm 
tjrbulence — i.e.,  turbulence  directly  associated  with  convective  or 
precipitating  cl^iuds  is  excluded.  Since  direct  turbulence  observations 
are  not  made  routinely,  o’ e  cannot  now  refer  to  an  archive  to  ootain 
data  for  making  a  climatological  study.  Instead,  an  alternate  method 
must  be  formulated  that  uses  various  types  of  fragnentary  information. 

The  general  plan  is  as  follows.  On  a  limited  number  of  occn.sions  In¬ 
strumented  airciaft  have  measured  turbulence  as  well  as  the  associated 
meteorological  conditions,  and  thus  provide  basic  irformit ion  concerning 
atmospheric  conditions  that  are  associated  with  turbuletce.  Next,  tur¬ 
bulence  frequencies  on  certain  days  were  computed  from  reports  of  tur¬ 
bulent  or  smooth  flight  collected  from  airline  pilots.  These  frequencies 
were  correlated  with  concurrent  meteorolcg ical  conditions  determined  as 
well  as  possible  from  standard  upper-air  data,  snd  regression  equations 
were  found.  The  final  step  was  to  apply  the  regression  equations  lietweer. 
turbulence  and  meteorological  factors  to  existing  climatological  recoros 
of  meteorological  conditions,  in  order  to  obtain  the  estimated  turbulence 
cl imatology . 

Numerous  difficulties,  some  anticipated  and  others  not  foreseen, 
are  inherent  in  this  indirect  approach.  Ttiese  will  be  dlsc’js»ed  in 
detail  later.  Conceptually,  the  simplest  wav  to  oty'a.n  a  tc-bulence 
climatology  would  be  to  equip  enough  aircraft  ti>  directly  determine  the 
turbulence  frequency  at  a  variety  of  places,  altitudes,  and  t.mes, 
however,  costs  would  be  prohibitive.  Another  approach,  which  we  have 
advocateo  but  which  is  still  untried,  would  be  to  develop  a  halloon- 
(x>rr,e  vurbuience  sensor  to  flv  as  part  of  the  standard  radiosonde  in- 
.strument.  vir  sepiiatcly.  If  an  inexpensive  sensor  proved  feasible. 

'he  turbulerue  ^ !  i  ma  t  <.i  1  ogv  up  to  the  altitude  li-.<t  of  balloon  flights 


(approximately  100,000  feet)  could  be  found,  and  correlations  could  be 
made  with  meteorological  factors  at  the  same  times  and  places.  A 
similar  Idea  Is  Involved  In  the  use  of  rising  uninstrumented  balloons 
tracked  by  hlgh-preclslon  radar  to  reveal  erratic  portions  of  the  se¬ 
quence  of  flight  coordinates  that  would  Indicate  turbulent  layers. 

This  subject  Is  discussed  In  Sec.  III.  Other  possible  methods  exist 
for  remotely  sensing  turbulence  using  el3ctromagnetlc  radiation,  but 
their  practicality  Is  not  established.  Such  methods  are  not  considered 
In  this  report. 

Among  designers  of  aircraft  and  missiles,  a  strong  Interest  In 
turbulence  ha:;  existed  for  years  (e.g, ,  Loving,  1966;  Houbolt,  1967). 

The  aeronautical  profession  has  perfected  complex  aircraft  instrumenta¬ 
tion  that  gives  information  on  turbulence  spectra,  probabilities  of 
gusts  greater  than  certain  threshold  values,  etc.  Using  this  approach 
to  determine  the  turbulence  environment  at  new  altitude  levels,  such  as 
the  stratosphere,  requires  that  an  aircraft  of  some  sort  be  used  as  the 
Initial  probe.  Generally,  the  number  of  flights  made  is  far  too  low  to 
indicate  v".riations  in  the  percentage  of  turbulent  air  to  be  expected 
under  different  weather  conditions,  places,  altitudes  and  seasons.  It 
Is  this  climatological  distribution  of  turbulence  that  is  of  primary 
interest  in  the  present  Investigation. 

This  report  does  not  summarize  or  review  the  voliuninous  literature 
on  turbulence  in  the  free  atmosphere;  however,  references  to  pertinent 
papers  are  given  at  appropriate  points.  In  Sec.  II,  the  different  data 
sources  that  we  used  are  described,  and  *heir  limitations  are  given. 
Section  III  is  entirely  devoted  to  a  description  of  our  attempt  to 
identify  turbulent  layers  from  paths  of  Rose  balloons  tracked  by  FPS-16 
radar.  Section  IV  contains  the  regression  equations  used  to  obtain  the 
estimated  turbulence  frequencies  of  Sec.  V.  (Examples  of  the  climatolo¬ 
gical  meteorological  data  used  are  given  in  the  Appendix.)  Recommenda¬ 
tions  aie  given  in  Sec.  VI,  A  scientific  report  concerning  the  first 
year's  work  was  issued  under  this  project  in  May  1966. 
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I I  DATA  SOURCES 
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As  discussed  In  Scientific  Report  1,  data  obtained  by  the  AFCRL 
B-47  research  aircraft  (which  measured  turbulence  intensity  objectively, 
and  also  recorded  concurrent  winds  and  temperature)  were  carefully  re¬ 
viewed.  In  recent  years  data  of  this  type  have  been  analyzed  by  a 
number  of  writers  including  Briggs  and  Roach  (1963),  Endlich  (1964), 
Endlich  and  McLean  (1965),  McLean  (1965),  and  Kao  and  Woods  (1966). 

From  such  information,  it  is  usually  possible  to  identify  detailed 
"mesoscale"  features  of  atmospheric  structure  that  accompany  turbulence. 
We  found  that  the  turbulent  gust  intensity  recorded  by  the  B-47  corre¬ 
lated  best  with  vertical  vector  wind  shear,  with  large  directional  wind 
shear,  and  with  deformation  in  the  horizontal  plane.  Correlations  be¬ 
tween  turbulent  gust  intensity  and  these  quantities  were  in  the  range 
from  0.5  to  0.8.  The  vertical  and  horizontal  structure  of  the  tempera¬ 
ture  field  appeared  to  be  less  important.  However,  the  mesoscale  me¬ 
teorological  patterns  that  produce  or  accompany  turbulence  are  not 
adequately  depicted  by  standard  upper-air  soundings.  Therefore,  rela¬ 
tionships  between  turbulence  and  meteorological  factors  are  invariably 
more  obscure  when  only  standard  meteorological  observations  are  available 
for  analysis.  The  correlations  of  this  report,  which  relate  standard 
meteorological  data  and  subjective  pilot  reports,  are  much  lower  than 
the  correlations  mentioned  above. 

InstruiTiented  gliders  and  aircraft  have  also  documented  the  severe 
turbulence  that  occurs  in  certain  portions  of  well  developed  mountain 
waves  (Kuettner,  1958;  Jones  and  Atnip,  1964).  Quite  recently,  insti’u- 
mented  U-2  aircraft  have  explored  stratospheric  turbulence  (Penn  and 
Pisinski,  1967;  Crooks,  et  al.,  1967).  These  data  apparently  show 
that  relationships  exist  between  turbulence  and  mesoscale  wind  features 
similar  to  those  found  earlier  in  the  troposphere. 

Since  the  variability  of  mesoscale  features  is  very  broad  while 
the  amount  of  hj gh-quality ,  detailed  aircraft  data  that  describe  them 
is  small,  it  is  difficult  to  generalize  from  the  aircraft  data. 
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Another  possible  source  of  Information  comes  from  balloons  tracked  by 
FPS-16  radar.  Flight  tracks  of  Jimsphere  balloons  had  been  obtained  by 
NASA  at  Cape  Kennedy^  Florida  (Scoggins,  1965),  and  of  Rose  balloons  by 
the  Air  Force  at  Vandenberg  AFB,  California,  for  use  In  measuring  accu¬ 
rate  winds  and  vertical  wind  shears  over  small  increments  of  height. 

It  was  thought  that  If  these  data  yielded  Information  on  turbulence, 
they  would  augment  other  turbulence  measurements,  particularly  In  the 
stratosphere.  From  the  analysis  of  ten  Jimsphere  flights,  the  approach 
appeared  feasible,  as  reported  by  Endllch  and  Davies  (1967).  As  a  re¬ 
sult,  a  major  effort  was  made  to  use  the  Rose  series  of  flights  from 
Vandenberg  for  turbulence  studies.  For  the  stratospheric  layers  of 
interest,  this  was  not  successful  for  reasons  discussed  further  in 
Sec.  III.  However,  the  basic  concept  of  detecting  turbulence  from 
erratic  portions  of  the  balloon  path  is  still  believed  to  be  useful, 
and  might  be  implemented  by  procedural  changes  described  in  Sec,  III, 

The  main  source  of  data  used  in  this  investigation  is  that  from 
airline  pilots,  who  have  reported  turbulence  severity  encountered  during 
commercial  flights  under  programs  organized  by  Colson  (1963,  1965,  1966) 
in  association  with  the  International  Civil  Aviation  Organization.  Each 
program  was  five  days  in  length.  Reporting  periods  were  held  in  February 
1963,  December  1964,  and  March,  June,  and  September  1965.  These  data 
have  the  major  advantage  that  they  are  numerous,  and  give  a  reasonably 
good  depiction  of  turbulence  within  the  airspace  over  the  United  States 
during  these  periods.  For  example,  in  the  March  1965  data,  there  were 
approximately  20,000  reports,  each  pertaining  to  a  100-mile  flight  seg¬ 
ment.  Turbulence  of  moderate  intensity  occurred  in  approximately  six 
percent  of  these,  while  severe  turbulence  occurred  only  in  0.3  percent 
of  the  cases.  The  average  wind  patterns  during  this  five-day  period 
were  by  no  means  typical  of  average  conditions  in  March,  as  can  be  seen 
from  upper-air  charts  for  this  interval  given  by  Endllch  and  Mancuso 
(1967).  The  same  is  true  of  other  five-day  periods.  Therefore,  one 
would  not  expect  the  pilot  reports  of  turbulence  for  such  brief  periods 
to  be  representative  of  seasonal  conditions — i.e.,  they  do  not  give  an 
approximate  turbulence  climatology.  But  by  matching  the  turbulence 


reports  to  concurrent  meteorological  conditions,  we  wished  to  obtain 
reliable  statistical  relationships  to  use  in  estimating  a  turbulence 
climatology. 

The  principal  disadvantage  of  the  pilot  reports  is  their  subjec¬ 
tivity  in  regard  to  the  severity  of  turbulence.  Another  difficulty  Is 
that  in  summer  considerable  turbulence  included  in  the  clear-air  category 
is  located  in  the  vicinity  of  convective  storms  or  downstream;  these 
cases  are  probably  an  indirect  result  of  the  convective  storms.  There¬ 
fore,  the  June  turbulence  has  different  relationships  to  meteorological 
factors  than  found  in  other  seasons. 

Colson  divided  the  airspace  into  volume  elements  approximately 
4000  feet  deep  and  bounded  by  2-1/2  degree  latitude-longitude  lines, 
and  determined  the  number  of  flights  during  twelve-hour  periods  that 
encountered  smooth  flight  and  turbulence  in  the  categories  light, 
moderate,  or  severe.  From  such  tabulations  furnished  to  us  by  Colson 
for  the  first  two  periods,  and  by  the  Federal  Aviation  Administration 
for  the  latter  three  periods,  we  computed  the  frequency  of  turbulence  in 
each  time  period  for  each  volume  element.  The  concurrent  meteorological 
conditions  (i.e.,  wind  speed  and  direction,  vertical  v.ctor  shear,  tem¬ 
perature,  lapse  rate,  horizontal  wind  shear,  deformation,  vorticity, 
divergence,  and  twelve-hour  vector  change  in  wind)  in  each  volume  were 
analyzed  by  computer  from  standard  rawlnsonde  data  at  synoptic  hours,  as 
described  by  Endlich  and  Mancuso  (1967).  Thus  a  complete  "turbulence 
report"  consists  of  the  turbulence  frequency  in  a  volume  element,  plus 
the  associated  wind  speed,  shear,  etc.  As  the  typical  flight  distance 
through  a  volume  element  is  approximately  lou  nautical  miles,  the  tur¬ 
bulence  frequencies  should  be  interpreted  as  the  likelihood  of  encoun¬ 
tering  turbulence  within  that  distance,  (This  distance  can  be  converted 
to  an  equivalent  flight  time  by  dividing  it  by  the  appropriate  aircraft 
groundspeed  , )  The  center  of  each  volume  element  is  represented  by  a 
grid  point.  Statistics  on  the  length  of  turbulent  patches  are  not  given 
by  these  data,  but  sources  such  as  Steiner  (1965)  and  Coy  (1967)  indicate 
that  in  approximately  50  percent  of  the  cases  the  length  is  less  than 
ten  miles. 
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To  use  the  pilot  reports  In  obtaining  regression  equations  for  the 
frequency  of  turbulence,  a  "climatology"  of  wind  speeds,  shears,  etc. 
at  each  grid  point  was  constructed  from  ten  values  of  each  factor  during 
each  five-day  period.  Then  a  multiple-regression  computer  program  de¬ 
veloped  by  M.  Gorfinkel  of  SRI  was  used  to  relate  the  frequency  (in  per¬ 
cent)  of  moderate  or  severe  turbulence  to  the  wind  speed,  vertical  shear, 
deformation,  lapse  rate,  etc.,  during  the  period.  Correlations  between 
turbulence  frequency  and  Individual  meteorological  factors  in  December 
and  March  are  given  in  Table  I.  Further  details  are  given  in  Sec.  IV. 

The  next  step  was  to  apply  the  regression  equations  to  existing 
climatological  records  of  pertinent  meteorological  factors  in  order  to 
estimate  the  turbulence  climatology  over  the  United  States.  It  was 
assumed  at  the  beginning  of  this  study  that  such  climatological  data 
would  be  available;  determination  of  new  climatological  descriptions  of 
the  upper  atmospheric  structure  was  not  contemplated.  Climatological 
data  concerning  wind  shear  over  the  United  States  were  given  by  Ratner 
(1958)  in  the  altitude  range  700  to  50  mb  (10,000  to  67,000  feet)  for 
the  period  1951-1956,  and  relied  upon  subjective  interpolation  to  com¬ 
pensate  for  missing  wind  observations.  The  amount  of  missing  data  was 
not  presented,  but  as  indicated  below,  must  have  been  large.  Further 
unpublished  wind  shear  statistics  mentioned  by  Crutcher  (1963)  had  been 
compiled  by  S.  Schamach  of  the  National  Weather  Records  Center  for  a 
ten-year  period.  These  were  generously  made  available  to  us  by  Dr, 
Crutcher.  They  are  for  levels  between  450  and  150  mb  for  the  period 
July  1948  through  June  1958,  The  statistics  are  comprehensive,  and 
obviously  excellently  computed  and  organized.  The  number  of  observations 
used  is  also  given;  no  subjective  interpolation  was  performed.  In  this 
compilation,  there  are  approximately  900  possible  observations  at  each 
station  for  each  season  (once  a  day  for  three  months  for  ten  years). 

The  actual  number  of  observations  shows  a  drastic  decrease  with 
altitude,  particularly  during  the  winter  season.  For  example,  at 
Washington,  D.C.  (one  of  the  better  stations),  532  wind  observations  In 
summer  reached  the  200-mb  level.  In  winter  at  200  mb,  the  corresponding 
number  Is  down  to  300.  .At  Ely,  Nevada,  the  numbers  of  wind  observations 
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at  200  mb  In  summer  and  winter  are  256  and  115,  respectively,  again  out 
of  900  possible.  Thus  the  samples  of  measured  winds  and  shears  are 
quite  small.  This  In  Itself  Is  not  of  overpowering  Importance;  the 
remaining  samples  might  be  representative.  But  from  the  nature  of  the 
observational  method,  we  know  that  wind  data  are  lost  most  frequently 
during  conditions  of  high  winds  when  elevation  angles  become  too  low 
to  permit  accurate  tracking.  Thus  the  remaining  wind  sample  Is  biased 
towards  low  speeds,  low  shears,  etc.  In  portraying  the  geographical 
and  latitudinal  variation  In  these  statistics  over  the  United  States, 
the  differing  numbers  of  observations  at  radiosonde  stations  enter 
Into  the  picture.  Although  we  have  looked  briefly  Into  the  possibility 
of  standardizing  the  statistics  on  the  basis  of  the  number  of  observa¬ 
tions  that  went  into  each,  no  simple  correcting  procedure  appears 
possible.  In  retrospect,  one  recalls  that  the  present  GMD-1  type  of 
wind  finding  equipment  was  introduced  at  some  stations  in  1954,  and  at 
most  others  in  196C.  Therefore,  the  number  of  observations  obtained  at 
each  station  evidently  depends  In  part  on  the  date  of  installation  of 
this  equipment. 

This  discussion  of  the  observations  that  go  into  available  statis¬ 
tics  of  ..pper  winds  and  shears  has  been  rather  lengthy  because  of  the 
importance  of  this  matter  to  the  present  problem.  The  reliability  of 
the  averages  and  standard  deviations  of  vertical  vector  shear  (used 
later  In  the  regression  equations)  is  difficult  to  judge,  and  may  be 
somewha*  questionable  for  the  present  purpose.  Typical  patterns  of 
the  magnitude  of  the  average  vertical  wind  shear  and  the  standard 
vector  deviation  of  shear,  as  read  at  grid  points  from  data  at  twenty- 
five  stations  given  by  Schamach,  are  shown  in  the  Appendix.  A  future 
remedy  appears  to  be  the  recomputat ton  of  such  statistics  using  recent 
data,  although  even  today  a  significant  number  of  winds  are  missed  in 
the  critically  important  Jet  stream  cases.  A  further  important  step 
that  can  be  taken  to  obtain  more  accurate  winds  and  shears  from  present 
equipment  Is  a  computer  computation  using  all  the  data  points  measured 
by  the  GMD-l  (instead  of  certain  points  one  or  two  minutes  apart),  as 
discussed  jy  Oanielsen  and  Duquet  (1987). 
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An  alternative  type  of  climatological  information  might  describe 
the  number  of  times  certain  subsynoptic  features  that  favor  t urbi.ilerice 
(ruch  as  Jet  streams  passing  through  sharp  troughs  and  ridges,  Jet 
1  routs  or  tropupauses  with  strong  wind  shear,  mountain  waves,  etc.) 
occur  as  a  function  of  geographical  location  and  time.  One  would  also 
need  to  know  the  probability  of  turbulence  associated  with  each  '"eteoro- 
logical  feature.  But  such  statistics  are  not  available,  and  anyway, 
these  phenomena  should  be  reflected  indirectly  in  means  and  standard 
ticv tat  ions  of  winds,  vertical  shears  and  temperature  lapse  rates. 

The  wind  stiear  data  of  Schamach  are  over  50-mb  intervals  i  jIow 
200  mb,  and  then  are  from  200  to  175  mb,  and  175  to  150  mb.  Thus  their 
avei-age  thicknesses  are  variable,  and  wc  wished  to  normalize  the  shear 
data  to  a  TOOO-fuot  thickness  as  used  in  analyzing  f.he  aircraft  data. 

A  well  known  teature  of  wind  shears  is  that  average  values  and  standard 
deviations  computed  for  thin  l„,ers  tend  to  be  larger  than  values  com- 
puicd  over  greater  distance  intervals  (D\oskin  and  Sissenwine,  1958). 
facto,  s  relating  shears  over  various  intervals  I'.  >'.e  been  given  by 
Essenwanger  (1963)  and  Ar  .endariz  and  Kider  (1966).  From  the  latter 
paper,  we  estimated  correction  factors  of  0.5,  0.67,  1.0,  and  1,3  for 
layers  1000,  2000,  TOGO,  and  6000  feet  in  thickness,  respectively,  and 
interpolated  for  other  iltickness  values.  For  e.xample.  an  average  shear 
computed  over  1000  feet  must  be  multiplied  by  0.5  to  normalize  it  to 
the  desired  -iOOO-toot  shear. 
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Ill  TURBULENCE  IDENTIFICATION  FROM  PATHS  OF  RISING  BALLOONS 


As  mentioned  in  the  Introduction,  quantitative  measure... v  nts  of 
turbulence  intensity  and  ai.sociated  vinds  and  shear  are  scarce.  There¬ 
fore,  we  investigated  the  applicability  of  data  from  Jimsphere  and  Rose 
balloons  tracked  by  the  highly  accurate  fPS-16  radar.  The  basic  concept 
is  simply  that  the  balloon  path  will  be  relatively  erratic  in  turbulent 
lavers  as  compared  to  laminar  portions  of  the  atmosphere  (Endlich  and 
Davies,  1£67).  A  preliminary  study  was  made  of  ten  Jimapti'^re  flights 
from  Cape  Kennedy,  and  these  appeared  to  indicate  the  existence  of 
several  turbulent  layers  under  seemingly  favorable  conditions.  With 
this  encouragement,  the  technique  was  applied  to  Rose  balloons  tracked 
by  Fre-16  radar  at  Vandenberg  Air  Force  Base,  California, 

As  discussed  in  the  refererced  paper,  data  processing  can  be  done 
in  terms  of  range  (R),  elevation  angle  (E),  and  azimuth  (A),  or  after 
converting  these  basic  measo'ements  to  coordinates  towards  east,  north, 
and  up  (x,  y,  z).  We  choose  to  use  the  former  three  coordinates.  Since 
the  variations  in  R,  E,  A  expected  due  to  turbulence  are  a  small,  high- 
frequency  component  superimposed  on  large  overall  trends,  for  turbulence 
to  be  detected  it  is  necessary  that  errors  due  to  radar  tracking  and  to 
self-induced  balloon  motions  be  relatively  small.  A  number  of  Rose 
balloon  flights  had  been  made  using  two  radars  located  side  by  .‘■Ide  to 
obtain  independent  measurements  of  the  balloon’s  co*.).~di  nates  versus 
time.  If  radar  tracking  errors  were  small,  plots  f ;om  the  dual  tracks 
should  be  essentially  identical,  ar  i  should  unequivocably  identify  tur¬ 
bulent  lavers. 

Computer  programs  were  written  to  utilize  -he  KPS-16  Rose  data. 
These  programs  are  as  fis)  lows  : 

(1)  The  first  program  converts  units  'f  R^  A  to  meters  and 

radians,  and  averages  five  ts.  1-second  ptunts  to  give  points 
at  0.5-second  intervals,  (T.'us  density  of  'nformation  is 
more  than  adequate  for  t.nc  present  purp.jsc,  as  the  balloon 
rises  at  a  rate  of  onlv  1  to  5  m  sec  )  The  time  sequence 


t  I 


of  points  is  then  examined  to  isolate  occasional  spurious 
points.  The  relevant  computation  is  made  between  adjacent 
points;  for  example,  if  Ir. 


R  . !  >  L  (where  L  is 

j-1 


approximately  50  m,  the  upper  limit  of  differences  in 


range  that  would  be  produced  by  a  lOU  m  sec 


wind 


in  0.5  seconds),  R  is  replaced  with  an  interpolated 
point.  Elevation  and  azimuth  are  treated  similarly, 

(2)  The  next  program  smooths  the  R,  E.  and  A  (separately) 

over  one-minute  intervals  using  a  triangular-shaped 
weighting  function.  Mean  wind  speed  and  direction, 
iteight,  and  balloon  ascent  rate  over  one-minute  inter¬ 
vals  are  computed  from  the  smoothed  values  (R^,  ^ ‘ 

Then  deviations  (R  -  R  ).  (E  -  E  ),  and  (A  -  A  )  are 

s  s  s 

computed  at  each  0.5  second,  and  a  plot  tape  is  written 
to  graph  the  deviations  as  a  function  of  time  using  a 
California  Computer  Company  plotter. 

(3)  Initially  it  was  also  planned  to  use  an  existing  computer 
program  to  compute  autocorrelations  and  spectra  of  the 
deviations,  but  this  was  not  done  for  reasons  given  below. 

Plotted  values  ot  (R  -  R  )  from  the  dual  tracked  flights  were 
graphed  as  shoh:.  in  Ftg.  1.  A  minute  in  time  corresponds  approxii,.ately 
to  1000  feet  in  altitude.  In  the  low  levels,  the  plots  from  the  two 
radars  generally  agreed  well,  and  indicate  self-induced  balloon  motions, 
and  occasional  turbulence  in  the  boundary  layer  of  the  atmosphere. 

Such  turbulence  is  shown  in  Fig.  1  for  the  interval  from  two  to  three 
m'nutes  after  launch.  oust  speeds  indicated  bv  the  slope  of  the  graphs 
are  1  to  2  m  sec  \  signifying  light  turbulence.  Other  similar  cases 
appart  tly  verify  the  basic  concept.  But  at  longer  ranges  and  low 
elevation  angles,  the  curves  the  two  radars  did  nut  generally  agree 
in  their  small-scale  features.  More  than  fifty  Rose  flights  were 
pri.)ces  sev'. ,  involving  a  fairli'  sizeable  effijrt  in  maniKiwer  and  computer 
costs.  In  this  sample,  tu-  clear-cut  cases  of  turbulence  in  ’.he  upper 
trojxisphere  i':  lower  s  t  ra  tt'.s  phere  could  be  identified  from  the  dual 
plot.s,  although  on-,  radar  or  the  I'ther  somet  iT4es  had  an  erratic  path 
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FIG.  1  GRAPHS  OF  RANGE  DEVIATIONS  MEASURED  BY  TWO  FPS-16 
RADARS  TRACKING  A  RISING  ROSE  BALLOON  IMMEDIATELY 
AFTER  LAUNCH.  Token  from  Test  6408  at  Vondenberg  AFB,  Californio, 
on  7  September  1966. 


FIG.  2  GRAPHS  range  DEV.aTiONS  MEASURED  BY  TwQ  FPS-16 
RADARS  ’■RACCSO  A  R  SING  ROSE  BALLOON  iN  THE  UPPER 
TROECSPHERc-  Ta^en  rrom  Test  6406  if  S  andenberg  AFB,  Califor’’  a 

on  T  September  1^6. 


suggesting  turbulence.  Close  agreement  was  not  found,  A  typical  plot 
(from  a  portion  of  the  same  flight  in  the  upper  troposphere)  is  shown 
in  Fig.  2.  The  duplication  of  the  two  curves  is  obviously  less  than 
in  Fig.  1. 

It  is  our  opinion  that  the  difficulties  are  due  to  tracking  uncer¬ 
tainties  that  increased  with  time,  particularly  at  long  ranges  where 
signal  strength  was  low.  These  uncertainties  could  not  be  removed  by 
the  editing  procedure  of  the  first  computer  program  mentioned  above. 
Also,  operational  modes  of  the  two  radars  may  not  have  been  optimum  or 
identical.  Thirily,  the  balloon  motion  may  be  partly  aerodynamic  in 
its  detailed  behavior.  Differences  between  the  motions  of  Rose  and 
Jlmsphere  balloons  have  been  noted  by  Scoggins  (1967), 

Unfortunately,  the  present  study  of  rising  Rose  balloons  tracked 
by  two  FPS-16  radars  was  not  successful  in  isolating  turbulent  layers. 
Still  we  do  not  consider  the  concept  discredited  since  the  Vandenberg 
flights  were  not  specifically  designed  to  obtain  measurements  of  tur¬ 
bulence.  We  suggest  a  trial  based  on  very  carefully  controlled  measure¬ 
ments  where  a  balloon  is  released  upstream  so  as  to  pass  over  the  two 
radars  at  a  favorable  (>  20  degree)  elevation  angle  and  a  relatively 
close  range.  The  assistance  of  a  search  radar  would  probably  be  needed 
to  aid  the  FPS-16's  in  locating  the  approaching  balloon.  At  Vandenberg 
AFB,  where  the  ocean  is  in  the  upstream  direction,  it  might  be  possible 
to  release  a  balloon  from  a  hel<copter.  In  spite  of  these  problems, 
further  experiments  with  FPS-16  tracking  of  rising  balloons  are  recom¬ 
mended  in  Sec.  VI. 
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IV  R2GRESSI0N  EQUATIONS  RELATING  TURBULENCE  FREQUENCY 
TO  METEOROLOGICAL  CONDITIONS 

A.  Winter,  Spring,  and  Fall 

As  mentioned  earlier,  by  far  the  most  numerous  turbulence  reports 
available  are  those  made  by  commercial  and  military  pilots  during 
special  five-day  reporting  periods  discussed  in  Sec.  II.  For  these 
periods  we  had  computed  turbulence  frequencies  In  volume  e.^ements 
covering  the  U.  S,  airspace  between  approximately  24,000  and  40,000  feet 
at  each  synoptic  hour  (00  and  12  GMT).  Also,  at  the  same  times  a  number 
of  meteorological  quantities  were  computed  from  objective  analyses  of 
standard  upper-air  observations  (see  Endllch  and  Mancuso,  19S7).  The 
selection  of  quantities  to  be  computed  was  made  on  the  basis  of  previous 
experience  with  data  from  many  sources,  and  on  our  estimate  of  what  can 
reasonably  be  determined  from  standard  rawinsondes.  [Opinions  on  this 
matter  differ;  for  example,  Moore  and  Krishnamurti  (1966)  advocate 
somewhat  different  quantities  than  ours.]  The  quantities  computed  for 
each  volume  element  included  zonal  and  meridional  wind  components  (or 
wind  speed  and  direction),  vertical  vector  shear  and  its  square,  vertical 
speed  shear,  temperature,  lapse  rate,  and  height.  From  the  winds,  vorti- 
city,  divergence,  deformation,  and  horizontal  shear  along  the  flow  and 
across  the  flow  were  computed.  The  values  of  all  these  quantities  are, 
of  course,  subject  to  various  errors  inherent  in  upper-air  data.  Fur¬ 
ther  derivatives  of  these  fields  (e.g.,  the  Laplacian  and  Jacobian) 
were  not  used  since  we  believe  such  terms  contain  an  unacceptable  level 
of  uncertainty. 

For  a  given  five-day  period,  a  "climatology"  was  constructed  for 
each  of  the  approximately  500  volume  elements  over  the  United  States. 

This  climatology  was  computed  from  the  ten  values  of  turbulence  fre¬ 
quency  and  ten  values  of  each  meteorological  quantity.  For  example, 
the  mean  vertical  vector  shear  and  standard  vector  deviation  of  shear 
were  computed  in  each  element.  Usually  about  50  volume  elements  had 
less  than  twenty  flights  during  the  period;  these  were  discarded.  Indi¬ 
vidual  correlations  of  turbulence  frequency  with  each  meteorological 
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quantity  were  determined  for  each  period  from  approximately  450  values 
of  each  variable.  Thus  the  sample  sizes  are  large.  Also,  a  multiple 
regression  program  (Gorflnkel,  1S67)  was  used  to  calculate  the  correla¬ 
tion  of  variables  with  each  other,  the  multiple  correlation,  and  the 
coefficients  of  linear  regression  equations.  This  program  also  gives 
the  standard  errors  In  the  estimate  of  each  coefficient  In  a  regression 
equation.  When  the  individual  correlations  are  all  relatively  lew 
(less  than  0.5),  as  In  these  data,  the  standard  errors  In  coefficients 
are  appreciable.  If  the  data  are  separated  into  sub-groups,  for  example, 
on  the  basis  of  altitude,  the  standard  errors  of  the  coefficients 
enable  one  to  Judge  whether  differences  in  the  coefficients  are  most 
likelv  real  or  random.  It  was  found  that  most  differences  among  sub¬ 
samples  (based  on  season  and  altitude)  were  random,  except  for  major 
differences  between  the  data  for  June  as  compared  to  the  other  periods. 

In  general,  the  turbulence  frequencies  (for  a  certain  period  at 
particular  points)  were  related  most  closely  to  vertical  vector  shear, 
and  somewhat  less  to  the  vertical  speed  shear.  The  probability  density 
function  of  vertical  vector  shear  is  not  known  precisely,  but  presumably 
approximates  a  circular  normal  distribution.  If  so,  its  properties  are 
well  described  by  the  vector  mean  and  standard  deviations.  Also,  it 
follows  that  the  mean  and  standard  deviation  of  the  speed  snear  could 
be  determined  from  the  vector  mean  and  standard  deviation  using  methods 
given  by  Brooks  and  Carruthers  (1954).  Both  types  of  statistics  (for 
vector  shear  and  speed  shear)  are  given  in  the  compilation  of  Schamach. 

Correlations  between  frequencies  of  moderate  or  severe  turbulence 
and  the  best  climatological  indicators  were  given  in  Table  I.  If  light 
turbulence  is  included  with  moderate  and  severe,  the  correlations  of 
Table  1  increase  by  amounts  of  approximately  0.1.  Thus  all  correlations 
obtained  from  these  data  are  considerably  lower  than  those  obtained 
using  data  from  special  research  aircraft,  as  described  in  Scientific 
Report  1.  Undoubtedly,  the  difficulties  discussed  earlier  have  ad¬ 
versely  affected  the  correlations.  Numerous  trials  were  made  using 
different  combinations  of  the  meteorological  quantities,  and  a  simple 
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linear  equation  based  on  the  mean  and  standard  deviation  of  vector  shear 
gave  a  multiple  correlation  as  high  as  when  additional  terms  were  used. 
The  relevant  empirical  equation  for  p(M),  the  percentage  frequency  of 
moderate  or  severe  turbulence  is 


p(M)  =  C(s)  +  2.0  X 


— 

— 

Magnitude  of 

Standard 

average  vertical 

-  0.5  y 

deviation  of 

vector  shear 

vector  shear 

»  — 

_  « 

(1) 


where  the  units  of  vertical  vector  shear  ir  the  brackets  are  10  sec  , 
and  C(s)  is  a  constant  depe  ^Ing  on  season,  taken  as  6  in  winter,  and 
5  in  spring  and  fall.  (These  values  are  in  accord  with  the  average 
turbulence  frequencies  for  samples  of  data  for  different  seasons.)  This 
equation  has  beon  applied  to  the  climatological  dr’^a  of  the  Appendix  to 
obtain  estimated  turbulence  frequencies  discussed  in  the  next  section 
of  this  report.  On  the  dependent  data,  the  multiple  correlation  is 
0.36.  The  standard  errors  in  the  coefficients  in  this  equation  (as 
given  by  the  program)  are  approximately  one-fifth  of  the  values  given. 

No  systematic  differences  among  these  coefficients  were  found  when  the 
turbulence  data  were  sub-divided  into  altitude  groups.  Also,  no  diurnal 
differences  in  turbulence  frequency  were  found. 

Since  climatological  wind  shear  data  are  available  only  for  the 
United  States  and  perhaps  a  few  other  industrialized  countries,  Eq.  (1) 
cannot  be  applied  in  other  geograpliical  areas.  However,  smoothed  maps 
of  mean  winds  and  the  standard  vector  deviation  of  winds  are  given  in 
various  sources  (such  as  Air  Force  Manual  SACM  105-2)  for  the  entire 
northern  hemisphere.  In  Table  I  we  saw  that  turbulen  frequency  is 
related,  although  weakly,  to  tliese  two  quantities.  Table  II  shows  them 
to  be  highly  correlated  with  each  other,  so  that  wind  speed  alone  will 
specify  turbulence  about  as  well  as  a  combination  oi  the  two.  The  pat¬ 
terns  of  strong  winds  over  the  United  States  given  in  SACM  105-2  compare 
in  broad  features  to  the  turbulence  frequency  estimated  from  Eq.  (1)  as 
one  would  expect  from  the  well  known  general  association  of  turbulence 
and  jet  streams.  The  rough  correspondence  of  turbulence  frequency  and 
average  wind  speed  can  be  seen  by  comparing  Fig.  5(c)  (see  Sec.  V)  with 
Fig.  A-l(c),  and  Fig.  5(d)  with  Fig.  A-l(f).  Actually,  the  regions 
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Table  II 


CORRELATIONS  BETWEEN  METEOROLOGICAL  FACTORS  BASED  ON 
OBJECTIVE  ANALYSES  FOR  SPECIAL  TURBULENCE  REPORTING 
PERIODS  IN  DECEMBER  1964  AND  MARCH  1965* 


Meteorological  Factor 

Factor 

1 

Factor 

2 

Factor 

3 

Factor 

4 

Factor 

5 

Factor 

6 

1. 

Average  vertical 
speed  shear 

1.0 

0.73 

0,54 

0.59 

2. 

Standard  deviation 
of  vertical  speed 
shear 

0,73 

1.0 

3. 

Magnitude  of  average 
vertical  vector  wind 

shear 

1.0 

0.  17 

0.40 

0.43 

4. 

Standard  deviation 

of  vector  wind 
shear 

0.17 

1.0 

0,24 

0.26 

5. 

Wind  speed 

0,54 

- 

0.40 

0.24 

1.0 

0.88 

6. 

Standard  deviation 

of 

0.59 

- 

0.43 

0.26 

0.88 

1.0 

*Each  correlation  is  based  on  approximately  900  pairs. 


estimated  to  have  the  highest  frequency  of  turbulence  lie  on  the  south 
side  of  the  centers  of  the  average  isotachs.  Nevertheless,  preliminary 
estimates  of  turbulence  from  average  wind  speeds  as  given  in  U.S.  Air 
Force  Manual  SACM  105-2  may  be  useful  for  certain  purposes.  Linear 
regression  equations  between  turbulence  frequency  and  average  wind 
speed  varied  considerably  among  subsamples  of  our  data.  Overall,  the 
regression  equation  that  relates  turbulence  frequency  to  average  wind 
speed  in  m  sec  ^  and  that  appears  to  give  reasonable  results  over  the 
United  States  is 

p(M)  =  C(a)  4-0.12  'i  [Average  wind  speed)  (2) 

where  C(a)  depends  on  altitude  and  is  taken  as  5  percent  at  400  mb, 

4  percent  at  300  mb,  and  3  percent  at  200  mb.  It  would  be  Interesting 
to  test  this  equation  in  other  geographical  areas,  particularly  Japan, 
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where  the  jet  stream  reaches  maxlhium  speeds.  However,  the  use  of  wind 
speeds  to  estimate  turbulence  frequency  is  less  •  .liable  than  the  use 
of  shear  statistics,  as  stated  previously. 

B.  Regression  Equations  for  Summer 

Compared  to  pilot  reports  for  several  wint:r  periods  that  are  dis¬ 
cussed  above,  reports  from  only  a  single  five-day  period  (June  9-14, 
1965)  were  available  to  represent  summer  conditions.  It  would  be  desir¬ 
able  to  verify  the  present  results  with  other  summer  data.  As  mentioned 
previously,  associations  of  turbulence  with  meteorological  quantities 
over  the  United  States  in  June  are  considerably  different  than  in 
winter.  The  mean  position  of  maximum  winds  in  summer  is  over  the 
Great  Lakes,  north  of  the  main  flight  routes.  The  vertical  wind  shear 
is  generally  much  weaker  than  in  winter.  i\^  a  large  degree  the  weather 
disturbances  over  the  United  States  are  convective  in  nature.  .Vluch  of 
the  convection  occurs  in  areas  where  vertical  shear  in  the  upper  tropo¬ 
sphere  is  small.  Even  the  turbuience  outside  convective  clouds  and  re¬ 
ported  by  pilots  as  clear-air  turbulence  may  have  been  indirectly  ol 
conve  tive  origin.  Probably  in  and  north  of  the  average  jet  stream  mi 
summer,  turbulent  frequencies  depend  primarily  on  wind  shear,  as  in  the 
other  seasons  over  the  United  States.  However,  the  June  data  show  nega¬ 
tive  correlations  of  turbulence  frequencies  with  average  vertical  shear, 
average  wind  speed,  and  the  standard  vector  deviation  of  wind.  Of  the 
quantities  computed,  only  relative  vorticity  has  a  small  positive  corre¬ 
lation,  as  shown  in  Table  III. 

Table  III 

CORRELATIONS  BETWEEN  PERCENTAGE  FREQUE.NCY  0.^  MODERATE  OR 
SEVERE  TURBULENCE  (FROM  PILOT  REPORTS)  .AND  CLIMATOLOGICAL 


Magnitude 
of  Average 


FACIXDRS,  FOR  JUNE  9-JJ,  1965 


Standard 

i'ecti.  i- 


Vertical  Vector  Deviation  of 
Wind  Shear  i  Wind  Shear 


-0.20  -0.15 


Each  correlation  is  based  on  approximately  400  pairs. 

The  best  regression  equation  fo'ind,  having  a  multiple  correlatio 
0.38,  is 


P(M) 


.  ».0  -f  1.7  X  y 


Magnitude  of  average 
vertical  vector  shear 

-1 

units  10  sec 


Standard  deviation 
of  vector  shear 

-3  -1 

units  10  sec 


'] 


-  0.5  X  Qrind  SpeedJ  4-0.1  v 

-1 


Standard  deviation 
of  wind  vector 


■] 


units  m  sec 


units  m  sec 


-1 


4-  0.03  (Relative  vortlclty) 

.  -1 
units  10  sec 


(3) 


Thl>7  equation  is  used  to  obtain  the  estimated  turbulence  frequencies 
for  summer  given  in  Sec.  V.  The  values  of  mean  wlna  speed  ind  direction 
in  summer,  and  the  standard  deviation  of  wind  were  taken  from  Air  Force 
Manual  SACM-105-2.  Vcrticlty  was  computed  from  the  mean  winds.  Typical 
examples  of  these  fields  are  shown  in  the  Appendix,  in  Fig.  A-3.  In 
these  summer  data  a  diurnal  effect  was  found,  such  that  computed  turbu¬ 
lence  frequencies  at  00  GMT  (evening)  should  be  Increased  by  15  percent, 
while  at  12  GMT  (morning)  the  probabilities  should  be  decreased  by  the 
same  amounts. 

C.  Mountain  Effects 

In  these  regression  equations,  no  separation  was  made  between 
mountainous  and  non-mountainous  terrain.  The  flight  frequency  is  such 
that  the  majority  of  data  pertain  to  non-mountainous  areas.  To  isolate 
terrain  effects,  grid  points  were  separated  into  mountainous  and  non- 
mountainous  categories.  Then  under  similar  meteorological  conditions, 
differences  in  turbulence  probability  between  the  two  categories  were 
found.  Inese  show  that  turbulence  frequencies  are  generally  higher  in 
mountainous  regions  than  would  be  indicated  by  a  regression  equation 
based  only  on  synoptic  parameters.  This  effect  can  be  adequately  rep¬ 
resented  by  increasing  the  turbulence  frequencies  in  mountainous  areas 
by  an  appropriate  amount.  The  mountainous  regions  of  the  United  States 
(generally  having  elevations  greater  than  1000  »)  are  shown  by  light 
shading  in  Fig.  3.  Within  these  art^as,  the  turbulence  frequencies 
given  in  Sec.  V  should  be  multiplied  by  a  factor  k^, 


taken  as  1.3  in 


FIG.  3  MAP  OF  PRINCIPAL  MOUNTAIN-WAVE  AREAS  OF  THE  UNITED  STATES. 

Within  th«  darker  shading,  turbulence  freguenci«r  of  subsequent  figures  should  be  multiplied 
by  kj  to  account  for  overage  mountain-wave  effects.  In  areas  of  lighter  shading 
the  multiplication  foctor  is  k^.  Values  of  k|  and  k^  are  given  in  Sec.  IV-C. 
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winter  end  1.15  in  spring  and  fall.  In  susuier,  no  such  correction  is 
needed.  It  is  known  that  the  largest  corrections  should  be  made  on  the 
lee  side  of  major  steep-sided  mountain  ranges  that  produce  prominent 
mountain  waves.  Such  areas,  shown  in  Fig.  3  by  dark  shading,  were  sub¬ 
jectively  estimated  from  topographic  maps  of  the  United  States,  and  in¬ 
clude  mountain  wave  regions  identified  by  Harrison  and  Sowa  (1966)  and 
Foltz  (1967),  To  account  for  wave  activity  in  these  areas,  turbulence 
frequencies  of  Sec.  V.  should  be  multiplied  by  k, ,  taken  as  1.6  in 

A 

winter  and  as  1.3  in  spring  and  fall. 

D.  Application  to  Other  Areas  and  Altitudes 

If  relevant  wind  shear  statistics  were  available  for  continental 
areas  cf  the  mid-latitudes  other  than  the  United  States,  it  would  be 
reasonable,  in  our  opinion,  to  apply  the  sane  regression  equations  In 
order  to  obtain  a  first  estimate  of  turbulence  frequency.  Over  the 
oceans  of  mid-latitudes  where  terrain  disturbances  are  lacking,  the 
regression  equations  probably  give  turbulence  frequencies  that  are  too 
large.  According  to  Clodman,  Morgan,  and  Ball  (1960),  turbulence  over 
oceans  is  an  order  of  magnitude  less  frequent  than  over  land.  However, 
recent  data  analyred  by  Colson  (1968)  Indicate  that  moderate  or  severe 
turbulence  over  oceans  is  about  0.7  times  as  frequent  as  over  land,  and 
that  the  same  sort  of  meteorological  conditions  (jet  stre'ims,  troughs, 
and  ridges)  are  important  for  both.  Thus,  a  simple  correction  of  this 
magnitude  might  be  used  in  applying  Eq.  (1)  to  the  oceans  of  mid¬ 
latitudes.  In  polar  regions,  one  might  expect  that  the  wintertime 
procedure  should  be  applied  in  all  seasons.  Perhaps  the  summer  regres¬ 
sion  equation  could  be  applied  In  the  tropics,  but  this  also  is  specula¬ 
tive  since  we  have  no  experience  with  aircraft  turbulence  data  for  tLese 
regions . 

it  must  be  borne  in  mind  also  that  the  most  comprehensive  turbulence 
data  (from  pilots)  have  been  obtained  in  ai.  altitude  range  restricted  to 
approximately  20,000  to  ^0,000  feet.  At  lower  altitudes,  the  regression 
equations  might  be  applicable  to  the  United  States  except  for  changes  in 
the  constant  values  to  accommcdate  greater  turbulence  amounts  due  to  the 
Influences  of  low-level  convection  and  terrain.  In  estimating  turbulence 
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Irequency  in  the  stratosphere,  one  is  faced  with  the  drastic  reduction 
in  available  information,  as  noted  by  Mitchell  (1966).  However,  the 
descriptions  of  U-2  measurements  of  turbulence  in  and  above  jet  streams 
Kiven  by  Penn  and  Pisinski  (1967)  indicate  that  the  wind  field  is  of 
predominant  importance,  with  a  possible  tendency  for  turbulence  to  be 
associated  with  intermediate  values  of  lapse  rate.  (The  contention 
that  stratospheric  turbulence  is  more  closely  related  to  temperature 
inversions  than  to  the  wind  field  appears  dubious  to  us.)  Due  to  the 
apparent  similarity  to  conditions  lower  down,  it  appears  reasonable  to 
apply  the  regression  equations  for  turbulence  to  the  lower  stratosphere. 
We  have  done  so  in  the  next  section,  subject  to  the  limitation  of  de¬ 
creasing  reliability  ol  standard  climatological  data  in  the  stratosphere. 
But  lor  altitudes  above  50,000  feet  a  new  range  of  problems  enters,  in 
that  future  aircraft  operating  at  these  levels  will  be  supersonic,  and 
will  be  affected  by  turbulent  eddies  of  longer  wavelengths  than  are  suo- 
soiuc  craft.  Spectrum  curves  of  turbulence  show  that  longer  wavelengths 
have  liigher  energy:  therefore  under  identical  atmospheric  conditions 
supersonic  planes  may  encounter  more  turbulence  than  would  conventional 
aircraft.  A  fairl\  large  amount  of  turbulence  encountered  in  a  small 
number  of  t lights  made  by  the  XB-70  aircraft  (Ehrenberger ,  196b)  ap¬ 
parently  supports  this  expectation.  For  this  reason  the  present  re¬ 
gression  equations,  if  applied  to  the  straivjsphere,  might  underestimate 
turbulence  probabilities  for  supersonic  aircraft.  This  problem  deserves 
further  investigation. 


V  AJi  ESTIMATED  TURBULENCE  CLIMATOLOGY  FOR  TOE  UNITED  STATES 


The  method  of  determining  the  turbulence  frequencies  presented  in 
this  section  was  described  in  Sec,  IV.  An  explanati'n  of  the  interpre¬ 
tation  and  use  of  the  fr  quencies  is  as  follows:  The  percentage  fre¬ 
quencies  p(M)  given  in  subsequent  figures  pertain  to  the  likelihood 
that  turbulence  of  moderate  or  severe  intensity  (as  evaluatec  by  pilots) 
will  be  encountered  during  a  100-mile  segment  of  a  flight  made  at  a 
particular  location  and  altitude.  The  average  value  of  p(M)  over  the 
United  Stats  is  shown  versus  altitude  in  Fig.  4.  If  light  turbulence 
is  combined  with  moderate  and  severe  to  give  p(T),  tl.e  overall  turbulence 
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frequency  (light,  moderate,  or  severe)  within  100-mile  sectors,  the 
pilot  reports  show  that  p(T)  =  3.5p(M)  in  winter  and  spring,  and 
p(T)  =  4.Sp(M)  in  summer  and  fall.  As  discussed  by  Endllch  and  Mancuso 
(1967),  the  frequency  of  encountering  the  same  intensity  of  turbulence 
in  a  flight  segment  shorter  than  100  miles  decreases  proportionately 
with  the  smaller  length.  Data  from  Steiner  (1965)  and  Endllch  (1965) 
indicate  that  the  frequency  of  encountering  moderate  or  severe  turbulence 
at  a  given  instant  is  less  than  two  percent.  This  amount  is  about  one- 
fifth  of  the  frequency  of  encountering  the  same  degree  of  turbulence  in 
some  portion  (or  possibly  all)  of  a  100-mile  segment.  Thus,  dividiLng 
the  frequencies  p(M)  by  five  will  give  p(m),  an  estimate  of  the  instanta¬ 
neous  frequency  of  encountering  moderate  or  severe  turbulence.  Similarly, 
dividing  p(T)  by  five  will  give  p(t),  the  instantaneous  likelihood  of 
encountering  light  or  greater  turbulence. 

In  summary,  the  turbulence  frequencies  that  follow  may  be  inter¬ 
preted  in  the  following  ways: 

(1)  The  frequencies  p(M)  of  this  section  pertain  to  the 
likelihood  of  encountering  moderate  or  se’ere  clear-air 
turbulence  within  a  100-mlle  flight  segment. 

(2)  To  obtain  p(T),  the  frequencies  of  encountering  light  or 
greater  turbulence  within  100-mile  flight  segments,  multiply 
p(M)  by  3.5  in  winter  and  spring,  and  by  4.5  in  suirjner 

and  fall. 

(3)  To  obtain  p(m),  the  likelihood  of  encountering  moderate  or 
severe  turbulence  at  a  given  instant,  divide  p(M)  by  5. 

(4)  Similarly,  to  obtain  p(t),  the  instantaneous  risk  of 
encountering  light  a  greater  turbulence,  divide  p(T) 
by  5. 

The  turbulence  frequencies  were  computed  at  5-degree  latitude- 
longitude  interrections  based  on  climatological  data  at  the  twenty-five 
stations  whose  locations  are  indicatf d  by  circles  in  Fig.  5.  Frequencies 
are  given  at  the  midpoint  of  those  altitude  layers  having  wind-shear 
data — l.e.,  at  425,  375,  325,  225,  187,  and  162  mb.  These  levels 
correspond  to  heights  of  approximately  22,000,  25,000,  28,000,  36,000, 


40,000  and  43,000  feet  respectively.  No  shear  data  were  available  at 
the  275-mb  (32,000-foot)  level.  Charts  for  winter  are  given  In  Figs. 
5(a)  through  5(f),  for  spring  In  Figs.  6(a)  through  6(f),  for  summer 
In  Figs.  7(a)  through  7(f),  and  for  fall  In  Figs.  8(a)  through  8(f). 

The  terrain  corrections  Indicated  by  Fig,  3  have  not  been  made  In 
Figs.  5-8;  this  must  be  done  by  the  reader  for  a  particular  locality. 
Similarly,  a  diurnal  correction  (Sec.  IV-B)  should  be  made  to  interpret 
Fig.  7  to  a  particular  time  of  day. 

These  turbulence  frequencies  are  believed  to  apply  to  problems 
where  turbulence  frequency  is  needed  by  season,  altitude,  and  place. 

They  provide  a  rough  basis  for  estimating  the  cverall  turbulence  expo¬ 
sure  associated  with  different  mission  profiles.  Operational  planning 
to  avoid  turbulence — for  example,  for  air  refuel: ng  missions--should 
obtain  some  guidance.  Also,  turbulence-seeking  p?  ^grams  may  be  directed 
to  areas  of  maximum  probable  exposure.  Other  applies 'ions  of  this  sort 
will  probably  be  found. 
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l?0*w  MO*  100*  90*  80*  7  0*w 


I20*W  110*  100*  90*  80*  70*  W 


FIG.  5(a)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  WINTER,  AT  425  mb 
(approx.  22,000  feet).  Average  percent  turbulence  =  7.0.  See  Sec.  V  for  an  interpretation 
of  these  frequencies. 


l?0*e  MO*  100*  90*  80*  7  0*W 


i20*e  MO*  100*  90*  80*  70*  e 


FIG.  5(b)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  WINTER,  AT  375  mb 
(approx.  25,000  feet).  Average  percent  turbulence  =  6.6. 
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i20*W  no*  100*  90*  80*  70*  W 


FIG.  5(c)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  WINTER,  AT  325  mb 
(approx.  28,000  feet).  Average  percent  turbulence  =  6.6. 


120**  no*  lOO*  90*  80*  7  0'w 


i20*VK  HO*  lOO*  90*  BO*  70*  W 


FIG.  5(d)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  WINTER,  AT  225  mb 
(approx.  36,(X)0  feet).  Average  percent  turbulence  =  .^.5. 


I20*w  Iio*  100*  90*  80*  70*  w 


FIG.  5(e)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  WINTER,  AT  187  mb 
(approx.  40,000  feet).  Average  percent  turbulence  =  3.8* 


120**  1 10*  100*  90*  80*  7  0*w 


i20*w  iiO*  lOD*  90*  80*  ’0*w 


FIG.  5(f)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  WINTER,  AT  162  mb 
(upprox.  43,000  feet).  Averoge  percent  turbulence  3.5. 


no*  100*  #0*  *0*  »0*W 


FIG.  6(a)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  SPRING.  AT  425  mb 
(approx.  22,000  feet).  Average  percent  turbulence  =  5.9.  See  Sec.  V  for  on  interpretation 
of  these  frequencies. 


i?0**  no*  lOO*  90*  #0*  ro*w 


FIG.  6(b)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  SPRING,  AT  375  mb 
(approx.  25,000  feet).  Average  percent  turbulence  =  5.5. 
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FIG.  6(c)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  SPRING,  AT  325  mb 
(approx,  28,000  feet).  Avera{]e  percent  turbulence  =  5.8. 


|?0*W  IIO*  100*  90*  SO*  7  0** 


FIG.  6(d)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  "TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  SPRING^  AT  225  mb 

(approx.  36, (XX)  feet).  Averoge  percent  turbulence  4.8. 
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FIG.  6(e)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  SPRING,  AT  187  mb 
(approx,  40,000  feet).  Averoge  percent  turbulence  ==  3.4. 
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FIG.  6(f)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUN  i  ckiNG  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  SPRING,  AT  162  mb 
(opprox.  43,(X)0  feet).  Average  percent  turbulence  3.0. 
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FIG.  7(q)  estimated  PROBABILITIES  (p«fc«nt)  OF  ENCOUNTERING  TURBULENCE  IN 
100-mile  sectors  over  the  united  STATES  IN  SUMMER,  AT  425  mb 
(approx.  22, OCX)  f«et).  Av«rag«  p«rc*nt  turb  'lanca  -  5.3.  S«e  S*c.  V  for  on  Intarpratotion 
of  th«so  frtquancio*. 
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FIG.  7(b)  ESTIMATED  PROBABILITIES  iporcorrt)  OF  ENCOUNTERING  TURBULENCE  IN 
lOO-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  SUMMER,  AT  375  mb 
(oppro.t.  25,000  f««r).  Avorog*  porcont  turbuUnc*  5.0. 
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FIG.  7(c)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCCJNTER'NG  TURBULENCE  IN 
100-MlLE  SECTORS  OVER  THE  UNITED  STATES  IN  SUMMER,  AT  325  mb 
(approx.  2B,(XX)  fo*f).  Avorog*  p«fC*r»t  turbuUnce  -  3.1. 


ilO*w  no*  '00*  »0‘  to*  rO** 


FIG.  7id)  EST'.MATED  PROBABILITIES  .p^rc.ofl  C-  ENCOUN'^ERING  TURBULENCE  IN 
'00-MiLE  SECTORS  OVER  THE  UNITED  STATES  IN  SUMMER,  AT  225  mb 
iopprox.  36,000  f***).  Av*ragt  p*rc*nt  furbui*nc*  3.3. 
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100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  SUMMER,  AT  187  mb 
(approx.  40,000  Avorap*  porcant  turbuionca  -  2.6. 


0*  'OC*  •C'  ’0*» 


<^!G.  7iT  ESTIMATED  PRCeASiLlTiES  parca-r?  OF  ENCOUNTER'NG  TURBULENCE  'N 
lOC-MILE  SECTORS  OVER  THE  UNITED  ST*TES  :N  SUMMER,  AT  162  mb 
30P»o«.  43.000  laat  .  Avirraga  parfart*  'u'bw'.arsea  2.2. 
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FIG.8(o)  ESTIMATED  PROBABILITIES  (ptretnt)  OF  ENCOUNTERING  TURBULENCE  IN 

lOO^ILE  SECTORS  OVER  THE  UNITED  STATES  IN  FALL,  AT  425  mb  (oppro*. 
22,000  ^Mt),  Avtrogt  ptretnt  turbultnct  =  5.3.  Stt  Stc.  V  for  on  inttrprttotion  of 
thtst  frtqutncits. 


HO’*  i  :•  lOo*  to*  to*  ’o*» 


FiG.  8ib^  ESTIMATED  PROBABILITIES  ptrct-fi  OF  ENCOUNTERING  TURBULENCE  IN 
lOO-MiLE  SECTORS  OVER  THE  UNITED  STATES  IN  FALL,  AT  3^5  -b 
jpores.  25,OOC  fttr  .  Avtrogt  ptrctn*  rur^.ltnet  5.3. 
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FIG.  8(c)  fcSTIMATk^D  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-M!Lt  SECTORS  OVER  THE  UNITED  STATES  IN  FALL,  AT  325  mb 
(a;-piox.  26,0'i0  fwt).  Average  psrcent  turSclence  =  5.6, 
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FIG.  8(d)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100  MILE  SECTORS  OVER  THE  UNITED  STATES  IN  FaLL,  AT  223  mb 
(approx.  36,00U  fees).  Average  percent  turbuier;ce  -  4,1. 
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FIG.  8(e)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  FALL,  AT  187  mb 
(approx.  40,000  feet).  Average  percent  turbulence  =  2.8. 
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FIG.  8(f)  ESTIMATED  PROBABILITIES  (percent)  OF  ENCOUNTERING  TURBULENCE  IN 
100-MILE  SECTORS  OVER  THE  UNITED  STATES  IN  FALL,  AT  162  mb 
(approx.  43,(X)0  feet).  Average  percent  turbulence  -  2.5. 
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VI  DISCUSSION  AND  RECOMMENDATIONS 


At  this  point  we  briefly  recapitulate  our  view  of  efforts  to 
estimate  the  turbulence  climatology  of  the  globe. 

Flights  by  instrumented  research  aircraft  provide  very  valuable 
basic  knowledge  about  turbulence  as  well  as  statistical  information 
for  design  purposes,  but  often  cannot  answer  important  questions  such 
as  the  following: 

(1)  What  is  the  three-dimensional  extent  of  a  turbulent  volume? 

(2)  What  are  the  special  conditions  in  the  generating  region 
of  the  turtulence? 

(3)  How  long  does  the  turbulent  area  last,  and  how  far  does 
it  travel  in  association  with  synoptic  meteorological 
features? 

(4)  If  the  mesoscale  meteorological  conditions  were  slightly 
different,  would  the  turbulence  still  have  been  present? 

Lack  of  answers  makes  generalization  from  limited  flights  quite 
difficult.  Future  flight  programs  should  attempt  to  obtain  inforr-cion 
concerning  these  matters  insofar  as  possible. 

Pilot  reports  have  been  of  major  importance  in  describing  turbu¬ 
lence  frequencies  at  altitudes  below  40,000  feet.  Since  equivalent 
information  is  not  available  in  the  stratosphere,  knowledge  of  strato¬ 
spheric  turbulence  will  tend  to  remain  incomplete  until  a  new  generation 
of  aircraft  has  been  built  and  flown.  At  that  point,  design  decisions 
will  have  already  been  made.  Mistakes  in  design  due  to  poor  knowledge 
of  turbulence  may  be  costly. 

Many  people  have  recognized  these  problems,  and  have  suggested 
that  remote  sensors  on  the  ground,  in  aircraft,  or  satellites  might  be 
able  to  detect  turbulence  or  associated  mesophenomona  (for  example, 
see  Reiter,  1967).  Probably  the  stratosphere  will  present  special 
difficulties  if  viewed  from  belc”.  Efforts  to  develop  remote  sensors 
are  being  carried  forward,  as  they  should  be.  However,  the  fact  that 
approximately  150  weather  balloons  rise  through  the  atmosphere  to  the 
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100,000-foot  level  In  the  United  States  every  day  has  been  ignored  in 
regard  to  the  possibility  of  their  carrying  direct  turbulence  sensors 
aloft.  The  sole  exception  in  the  United  States  is  the  very  interesting 
work  reported  by  Anderson  (1957)  using  the  gustsonde;  since  that  time 
technology  has  improved  greatly.  Certainly  on  low-level  towers,  direct 
measurements  of  turbulence  predominate  over  indirect  sensing,  and  this 
might  be  taken  as  a  partial  guide  for  the  stratosphere.  Instrumental 
problems  related  to  balloon-borne  sensors  are  unknown,  but  might  be 
much  simpler  to  solve  than  the  similarly  uncertain  problems  of  remote 
sensors.  Development  of  a  balloon-borne  turbulence  sensor  should  be 
undertaken  immediately  as  a  possible  means  of  augmenting  other  types 
of  data.  Also,  such  a  sensor  might  be  useful  prior  tc  flights  oi  in¬ 
strumented  aircraft  in  locating  favorable  atmospheric  layers  to  be 
probed  for  turbulence. 

As  discussed  in  Sec.  Ill,  special  trials  are  recommended  to  deter¬ 
mine  whether  the  FPS-16  Rose  system  can  detect  turbulence  in  the  strato¬ 
sphere.  Of  course,  this  can  only  be  done  at  those  few  sites  having 
FPS-16  radars. 

Since  the  presently  available  statistics  of  wind  shear  over  the 
United  States  are  based  on  data  obtained  prior  to  1960  that  are  unrep¬ 
resentative  of  jet-stream  conditions,  it  would  be  desirable  to  recompute 
them.  This  is  a  fairly  large  but  straightforward  task  that  could  be 
done  from  recent  higU-altitude  wind  measurements  made  by  OlD-l  (or 
equivalent)  equipment.  Appreciable  differences  between  the  new  wind- 
shear  climatology  and  earlier  values  would  dictate  that  the  estimated 
turbulence  frequencies  of  Sec.  V  be  redetermined. 


APPENDIX 


SELECTED  CLIMATOLOGICAL  DATA 
FOR  THE  UNITED  STATES 
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I20*W  no*  (00*  90*  90* 


(a)  Magnitude  of  the  mean  shear  vector  at  325  mb  in  winter,  units  10  sec 

I20*v»  1(0*  100*  90*  BO*  rO*W 


(b)  Standard  deviation  of  the  shear  vector  at  325  mb  in  winter,  units  10  ^  sec  ’ 

FIG.  A-1  STATISTICS  ON  THE  VERTICAL  WIND  SHEAR  VECTOR  (normoliied  too 
4000-foot  thickness)  IN  WINTER,  BASED  ON  SCHAMACH'S  DATA 
FOR  STATIONS  INDICATED  BY  SMALL  CIRCLES,  AND  WIND  SPEEDS  FROM 
U  S.  AIR  FORCE  MANUAL  SACM  105-2 
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(f)  M«on  wind  «p«*d  of  225  mb  in  wintvr,  unit*  m  tvc 


FIG.  A-1  (Continued} 
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(q)  Moqnttudtf  of  the  mean  sheer  vector  ot  325  mb  ir>  spring,  units  10  sec 

0*  •^-*  9C*  SC*  ’0* 


ib>  Stondnrd  deviOtion  of  the  l^eor  vector  ot  325  mb  m  tprtr»g,  un.*i  1C~  '  tec* 

FiG  A-2  STATI:>T1CS  ON  THE  VERTICAL  WIND  SHEAR  VECTOR  norwo  I.  t*d  *0  0 
4000-foor  IN  SPRING  BASED  0  4  SC-^AMACH'S  data 

FOR  STATiONS  INDICATED  BY  SMALL  CIRCLES 
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\-m)  Mo9nitu<W  of  moon  ifioor  voctor  at  187  "»b  m  spfn^,  units  10  '  sac 

JO**  :•  ic^ 
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*'  Stanford  ^»‘ation  of  ?ba  sha-wr  *ac^or  or  18!*  -rt  10  sac 

"  *0  A  O  Co^r  vaod 
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(a)  Magnitud*  of  the  moon  shoar  voctor  ot  325  mb  in  summer,  units  10~^  sec"' 
i20‘e  no*  too*  90*  eo*  ro*w 


FIG.  A-3  STATISTICS  ON  THE  VERTICAL  WIND  SHEAR  VECTOR  (normalized  to  a 
4000-foot  thickness)  IN  SUMMER,  BASED  ON  SCHAMACH’S  DATA 
FOR  STATIONS  INDICATED  BY  SMALL  CIRCLES,  AND  WIND  SPEEDS 
FROM  U.S.  AIR  FORCE  MANUAL  SACM  105-2 
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(c)  Mean  wind  speed  at  325  mb  in  summer,  units  m  sec” 

I20*W  no*  100*  90*  80* 


I20*W  no*  lOO*  90*  80* 


(d)  Stondord  deviotion  of  the  wind  vector  of  325  mb  in  summer,  units  m  sec 

FIG.  A-3  (Continued) 


55 


(•)  Relative  vorticity  computed  from  mean  winds  at  325  mb  in  summer,  units  10  ^  sec 
I20*W  no*  '00*  90*  80‘  rO‘W 


(f)  Magnitude  of  the  meon  shear  vector  at  225  mb  in  summer,  units  10  ^  sec 

FIG.  A-3  (Continued) 
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(g)  Standard  deviation  of  the  shear  vector  at  225  mb  in  summer,  units  10  ^  sec 

i2o*w  iio*  (00*  $0*  eo*  ^o*w 


(h)  Magnitude  of  the  mean  sheor  vector  ot  187  mb  in  summer,  unitt  10  tec 

FIG.  A-3  (Continued) 
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(a)  Magnitude  of  the  mean  shear  vector  at  325  mb  in  foil,  units  10”^  $ec“' 

HO*  lOO*  90*  eo*  7 


(b)  Standard  deviation  of  the  shear  vector  ot  325  mb  in  foil,  units  sec 

FIG.  A-4  STATISTICS  ON  THE  VERTICAL  WIND  SHEAR  VECTOR  (normoliied  to  o 
4000-foot  thickness)  IN  FALL,  BASED  ON  SCHAMACH'S  DATA 
FOR  STATIONS  INDICATED  BY  SMALL  CIRCLES 
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(e)  Magnitude  of  the  mean  sheor  vector  of  187  mb  in  fall,  units  10'^ 

i20*»  "0*  lOO*  90*  80* 


f)  Stondard  deviation  of  the  sheor  vector  ot  187  mb  m  foil,  units  10  ^ 


F!G.  A-4  Concluded' 
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